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GRAPHICS  DISCLAIMER 

All  figures,  graphics,  tables,  equations,  etc.  merged  Into  this 
translation  were  extracted  from  the  best  quality  copy  available. 


A  STUDY  OF  THE  ROTOR  WAKE  IN  NAP-OF- THE- EARTH 


He  Chengjian  and  5ao  Zheng 


(Nanjing  Aeronatuical  Institute) 


SUMMARY  This  article  does  research  into  the  primary  aerodynamic 
problem  with  *ap-of-the-earth  flying-the  ground  vortex  phenomenon. 
Through  an  analysis  of  rotor  wake  close  to  the  ground,  we  have  set  up 
models,  formulas,  and  methods  for  the  calculation  of  its  height  off 
the  ground,  its  longitudinal  position,  and  its  strength.  Moreover,  we 
have  used  these  to  calculate  an  induced  speed  distribution  at  the 
rotor  disc.  When  we  compared  the  results  of  calculations  and 
experimental  values,  we  arrived  at  relatively  good  agreement. 
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I.  Introduction 


We  had  done  research  early  on  concerning  the  ground  effects  of 
helicopter  rotors.  In  the  past,  nap-of-the-earth  flying  for 
helicopters  had  been  confined  only  to  the  short  term  phase  of  taking 
off  and  landing.  Research  on  ground  effects  had  been  primarily 
concentrated  on  the  aspects  of  rotor  tensile  strength  and  the  power 
required. 

In  recent  years,  due  to  the  requirement  for  the  use  of 
nap-of-the-earth  flying  by  military  helicopters  in  combat  zones,  there 
has  been  renewed  interest  in  rotor  ground  effect  problems  in 
helicopter  technology  circles  fl] .  Continuous  flight  close  to  the 
earth  at  low  speeds  causes  helicopters  to  be  placed  in  a  special 
aerodynamic  environment  and  produces  a  series  of  new  problems.  Among 
these  one  finds  involved,  during  periods  of  nap-of-the-earth  flight, 
helicopter  equilibrium,  handling  characteristics,  flight  capabilities, 
and  interference  from  surrounding  flow  fields,  as  well  as  a  number  of 
other  similar  helicopter  aerodynamic  problems. 

Rotor  wake  research  is  the  foundation  of  the  aerodynamics  of 
helicopter  nap-of-the-earth  flight.  It  is  only  when  one  can  set  up  a 
rotor  wake  model  which  is  objective  and  realistic  that  it  is  then 
possible  to  realistically  predict  the  aerodynamic  characteristics  of 
nap-of-the-earth  flight. 

In  classical  rotor  ground  effect  theoretical  research  12],  the 
rotor  wake  is  taken  as  fixed  in  order  to  describe  the  vortical  system 
model.  Making  use  of  the  projected  image  method  to  solve  for  the 
ground  effect  induction  coefficient,  it  is  possible  to  make 
reasonably  good  predictions  of  rotor  capabilities  in  hover  and  during 
forward  flight  at  relatively  high  speeds.  However,  application  of 
this  technique  is  limited  in  the  low  speed  ranges  encountered  during 
nap-of-the-earth  flight.  The  reason  for  this  is  that  no  account  is 
taken  of  the  relatively  large  wake  deformation  and  surface  vortical 
phenomena  at  such  times.  During  the  last  few  decades,  people  have 
achieved  a  very  rich  research  understanding  of  the  rotor  aerodynamics 
of  nap-of-the-earth  flight.  This  is  particularly  true  of  the  research 
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work  done  by  the  Boeing-Ford  Helicopter  Company  (5)  and  Princeton 
Univerity  (4,5).  This  work  makes  a  systematic  discussion  of  ground 
vortical  phenomena,  and,  from  experimentation,  presents  their 
processes  of  occurrence,  development,  and  disappearance.  Princeton 
University  also  measured  and  recorded  the  rules  governing  changes  in 
the  locations  of  rotor  flow  fields  and  ground  vortices  during 
nap-of-the-earth  flight.  This  is  precious  information  for  wake 
research.  Reference  [5]  carried  out  theoretical  research  into  rotor 
wakes,  ground  vortices,  and  mutual  interference  with  the  ground 
surface.  It  achieved  relatively  good  results.  Their  work  delved 
deeply  into  research  concerning  rotor  aerodynamics  during 
nap-of-the-earth  flight  and  laid  a  valuable  foundation. 


In  order  to  be  able  to  carry  out  predictions  of  the  aerodynamic 
characteristics  of  rotors  during  nap-of-the-earth  flight,  it  is 
necessary  to  solve  a  key  problem  which  is  the  calculation  of  the 
position  and  strength  of  ground  vortices.  The  research  which  previous 
people  have  done  has  still  not  provided  a  method  which  is  practically 
feasible  and  theoretically  complete.  This  article  presents  a  method 
for  calculating  the  strength  and  position  of  ground  vortices  and  rotor 
wakes  as  well  as  induced  velocities.  Because  of  this  method,  the 
calculated  results  and  the  experimentally  measured  data  from  Princeton 
University  (5)  mutually  agree. 


II.  Surface  Vortices 

In  recent  decades,  research  from  outside  China  has  pointed  out 
the  phenomenon  of  ground  vortices  which  appears  particularly  in 
nap-of-the-earth-flight.  Experimentation  [4]  discovered  that,  from  the 
appearance  to  the  disappearance  of  ground  vortices,  the  corresi onding 
rate  of  rotor  advance  lay  in  a  very  narrow  range  from  ^—0.02  to 
0.05.  At  this  time,  the  induced  velocity  distribution  at  the  rotor 
blade  disc  changes  abruptly  (Fig.  l).  It  creates  a  deformation  in  the 
rotor  aerodynamic  characteristics.  Its  pitching  moment  follows 
changes  in  |i  as  shown  in  Fig.  2. 
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Pig.  1  The  Relationship  Between  the  Axial  Component  of  the  Rotor  Dis 
Induced  Velocity  Kj  and  Its  Corresponding  Radial  Position  ».r/* 

[5] 

Pig.  2  Curve  Describing  Changes  in  the  Pitching  Moment  In  Terms  of 
the  Rate  of  Advance  M  17]  (  •.  nate  of  Slope  of  the  Lift  Curve,  '«  Degre 
of  Rotor  Reality)  3-  Influence  of  the  Ground  Vortex 

The  ground  vortices  are  an  accumulation,  in  the  vicinity  of  the 
ground  surface,  of  the  rotor  wake  vortex  (principally  the  rotor  tip 
vortex)  under  the  triple  mutual  influences  of  downward  rotor  wash 
flow,  the  oncoming  flow  encountering  the  flight  and  the  ground 
surface.  The  rotor  wake,  due  to  the  blocking  effect  of  the  ground 
surface,  reduces  speed  and  finally  turns  in  the  direction  of  the  flow 
It  is  dispersed  along  the  ground  surface  in  all  directions.  If  the 
speed  of  forward  flight  is  very  low,  the  forward  section  of  the  rotor 


wake  flow  still  has  a  relatively  large  horizontal  forward  velocity 
along  the  ground.  Because  of  the  continuity  of  the  flow  movement,  the 
wake  flow  disperses  and  reduces  speed.  At  a  certain  point  in  the 
forward  area,  it  mates  with  the  corresponding  wind.  The  supported 
wake  vortices,  then,  accumulate  at  this  spot  forming  the  ground 
vortex.  This  vortex  has  an  obvious  effect  on  the  surrounding  flow 
fields  and  the  aerodynamic  characteristics  of  rotors. 

1 .  Ground  Vortex  Vertical  Position  C, 

Since  ground  vortices  are  the  result  of  the  mutual  effects  of  the 
three  factors  of  rotor  wake  flow,  oncoming  flight  flows  and  the  ground 
surface,  the  position  of  ground  vortices  can  be  solved  on  the  ;he 
basis  of  the  principle  of  the  superposition  of  movements. 

After  rotor  wake  flow  collides  with  the  ground,  it  assumes  the 
forward  wake  flow  movement  of  the  wake  vortex  and  reverses  direction 
in  terms  of  the  oncoming  flight  flow.  The  horizontal  velocity  of  gas 
flow  close  to  the  ground  is  taken  from  the  superposition  of  the  two. 

At  the  point  where  their  flow  speeds  are  the  same  (their  combined 
speed  is  zero)  ,  the  wake  flow  turns  over  and  upward  forming  the  ground 
vortex.  The  vertical  position  of  this  vortex  can  be  solved  for  using 
the  classical  vortex  system. 

If  we  take  the  rotor  coordinate  system  as  our  reference  system 
(Fig.  3),  the  vertical  speed  of  the  surface  flow  is 


K,—  V  +  vH 


(D 


In  this  equation,  V  is  speed  of  flight.  vt,  is  the  rotor  wake 
vortex  and  its  horizontal  induced  velocity  as  projected  on  the 
vertical  plane  of  symmetry  rising  from  the  ground  surface.  On  the 
basis  of  classical  rotor  vortex  flow  theory  [6]  and  projection 
methods,  one  has 


(2) 


„  __  R 

f  * «  dr  r 

CDS0(#,COSX  —  A) 

_  _coeJcosX_'j 

m  2* 

J  «  rfl'l 

Rp(  Rp+Bp) 

Rn+Rfi  J 

In  this  equation  ~jf  is  the  amount  of  circulation  for  a  unit  length 
along  the  direction  of  the  center  line  of  the  vortex 

Rp~  V R* + tOOs6 
Bt—  (£? ,+ /?cos0)sinX — AoobX 

t/Rx  +  ( G»+ AtgX  )* + 2R(Gt + AtgX)cos9 
Btt  ■»  (  G . + RcosO  +  At  gX  )sin  X 


At  the  place  where  the  ground  vortex  is  formed  one  should  have  V.—  0 , 
that  is, 

V — (3) 


This  equation  determines  the  vertical  position  G«  of  the  ground 
vortex.  At  this  point,  it  is  necessary  to  point  out  that  the  wakes  in 
the  various  equations  above  must  be  used  with  correction  values  for 
the  angle  of  incline  X  of  the  vortical  core. 


Pig.  3  Rotor  Coordinate  System  1 .  Ground  Vortex 


Fig.  4  Curve  Expressing  Changes  in  the  Vertical  Position  c.  of  the 
Ground  Vortex  as  a  Function  of  »*  (  is  the  Overall  Distance) 


G,  is  an  important  parameter  describing  ground  vortices.  Fig.  4 
presents  changes  in  G, as  a  function  of  the  rate  of  forward  rotor 
movement.  In  all  situations,  <?.,  is  always  a  monotonic  reduction 
function  of  Mj.  The  diagrams  show  the  effect  of  the  overall  distance  "0,. 
on  the  position  of  the  ground  vortex.  The  larger  the  overall 
distance  is,  the  greater  are  the  momentum  and  the  kinetic  energy  of 
the  downward  wash.  The  ground  vortex  then  appears  in  a  place  even 
farther  to  the  front  of  the  rotor. 

Although,  when  we  solved  for  the  vertical  position  G,  of  the 
ground  vortex,  we  did  not  consider  the  influence  of  the  ground  vortex 
itself,  it  is  true,  however,  that  classical  vortical  systems  already 
calculate  in  that  portion  of  the  near  ground  wake  which  goes  to  form 
the  ground  vortex.  Because  of  this,  we  get  a  certain  compensating 
factor.  From  a  comparison  of  the  calculated  values  presented  in  Fig. 

5  and  experimentally  measured  values,  we  can  see  that  the  degree  of 
agreement  between  the  two  is  quite  satisfactory. 
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Fig.  5  A  Comparison  Between  Values  Calculated  for  G*  on  the  Basis  of 
the  Methods  of  This  Article  and  Experimental  Values  1  .  Experimental 
Values  2.  Calculated  Values 

<«)  =  10*,  (  ft  )•».»»  =  12*, 

2.  The  Height  Off  the  Ground  of  Ground  Vortices 

According  to  what  we  see  in  Fig.  6,  the  rotor  wake  flow  flows 
down  at  an  angle  toward  the  ground  at  velocity  V,  .  In  order  to 
simplify  this,  let  us  assume  that,  after  the  wake  flow  contacts  the 
ground,  it  divides  into  two  longitudinal  flow  movements  along  the 
ground.  The  thickness  of  the  forward  branch  is  .  The  thickness 

of  the  rearward  branch  is  *»•  .  When  we  pay  attention  to  applying  the 
principles  of  momentum  and  continuity  equations,  we  make  a  three 
dimensional  correction  to  the  wake  flow.  If  we  do  not  figure  in  the 
ground  friction  force,  it  is  possible  to  write  an  '*u  direction 
momentum  equation  for  the  control  body  shown  in  trie  diagram. 


Pig.  6  Derivation  of  Momentum  Equation  for  Control  Body  1 
Body 


Pi~P.  _ 
PQ  XRX  Gt 


VisinX 


The  continuity  equation  is 


.  Control 


(4) 


<7«y  +  fc,v, 


(5) 


If  we  make  equations  (4)  and  (5)  simultaneous,  and,  we  take  into 
consideration  the  energy  losses  which  occur  when  the  rotor  wake 
collides  with  the  ground  surface,  we  can  go  through  mathematical 
operations  and  obtain 


<?*- 


sinX) 


Rt  [~f~ (  ^!  +  8in*X,)  +  /C'lsinX  J 


3  29 


(6) 


In  this 


x 
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B * 


We  should  pay  attention  to  the  fact  that  the  rotor  blade  tip 
vortex  is  only  distributed  on  the  wake  flow  surface.  Because  of  tnis, 
the  thickness  of  the  branch  flow  in  the  forward  direction  after  the 
wake  flow  contacts  the  ground  here  is  then  the  height  of  the  ground 
vortex  which  we  solve  for.  Fig.  7  shows  a  comparison  of  calculated 
values  and  experimental  values.  The  two  agree. 


Fig.  7  A  Comparison  of  Experimental  Values  and  Values  for  ic*  ns 
Calculated  According  to  the  Method  of  This  Article  1 .  Experimental 
Values  2.  Calculated  Values 
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Let  us  assume  that  the  wake  vortex  which  avoids  the  rotcr  t-.ade 
tip  takes  the  form  of  a  vortical  ring  and  moves  downward  following  tn 
wake  flow.  At  a  place  a  certain  distance  above  tr.e  ground,  it  change 
and  is  only  dispersed  in  a  horizontal  plane.  Tne  forward  section  of 
the  ro*  'r  tip  vortical  ring  moves  slowly  due  to  its  windward  position 
The  rear  section  moves  a ",  high  speed  because  it  is  moving  with  the 
wind.  In  oruer  to  simplify  calculations,  take  the  wake  which  is 
s.nown  close  to  the  ground  by  the  planar  ring-shaped  vortex  of  the 
continuity  distribution  for  the  amount  of  circulation  v  Fig.  6 )  and 
assume  tna  +  its  circulation  distribution  is 

SO  (P</?) 

(  1  +  P OOB0j  ]  (  /?<  P  </?,) 


In  t  n  i  .■  h-VJ\\ 

Tnis  type  of  distribution  reflects  the  principle  of  longitudinal 
distriburion  of  the  amount  of  circulation  in  close  waxes  which  says 
tnat  the  front  is  large  and  the  rear  is  small. 

In  Fig.  8,  this  type  of  planar  vortex  system,  due  to  trie  mutual 
effects  of  vortex  rings,  has  edges  wnicn  are  unstatu-e,  ana  viu  roi. 
up.  The  forward  peak  of  tnis  vortex  system  is  at  tne  place  where  the 
oncoming  flow  of  flight  and  the  rotor  wake  come  together.  Under 
limits  imposed  by  tne  ground  boundary,  tne  waxe  flow  accumulates  at 
this  point,  turns  upward  and  rolls  to  form  a  ground  vortex. 

The  ground  vortex,  on  the  one  hand,  receives  a  continuous 
replenishment  from  the  rotor  tip  vortex  in  tne  rotor  down  f.ow.  n 
the  other  hand,  it  dissipates  its  amount  of  vortical  ity  tc  tne 
8urroundings  due  to  gas  flow  viscosity.  tonen  the  processes  of  vertex 
building  and  dissipation  reach  a  dynamic  e  }ui  i  or  i  urn ,  trier.,  tr.«- 
surface  vortex  is  set  up  and  maintained  at  a  certain  strength. 


In  this  9*tt  are  respectively  the  components  of 

induced  velocity  -.long  the  ’-x.axis  caused  at  tnat  point  in  the  rotor 
wake  in  response  to  changes  in  the  vortical  system  and  its  projection 
the  planar  vortical  system  pro j e c - i on ,  and  the  ground  vortex 
projection.  Since  tn«r  ground  vortex  and  th*-  rotor  are  stationary 
relative  to  each  other,  then.  K*—  o  .  From  this,  it  is  possicic 
to  obtain  the  strength  coefficient  for  tne  ground  vortex,  whicn  is. 


P,-Fl+tGl 

+  (k— (7,)tgX 


r.~ 
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for  'rr«T»'«  in  Terns  o'  Charges  in  • 
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Fig.  9  presents  the  principles  governing  the  average  values  of 
the  ratio  between  surface  vortex  strength  an  oln.de  tip  vortex 
strength  in  response  to  changes  in  u-  .  It  also  presents  tne 
influences  of  the  heignt  cf  the  rotor  off  the  ground  and  its  overall 
distance.  One  can  9ee  that  pro'jn  vortices  exist  ormy  in  a  very 
narrow  range  of  values  of  »  .  Moreover,  when  m  a as  a  certain  value, 
there  is  a  maximum  value.  Tne  greater  h  is  or  tne  smaller  0».?»  is, 
the  smaller  is  the  peax  value  for  '  (rVr,),  .  moreover,  values  of  ** 

which  correspond  to  these  peak  values  are  also  small.  These 
principles  of  change  and  the  patterns  of  experimental  observations 
agree  with  each  other.  It  would  be  reasonable  to  make  an  analysis  of 
this.  Unfortunately,  there  are  still  no  data  from  experimental 
measurements  with  which  to  make  comparisons. 

According  to  the  position  and  strength  of  ground  vortices  as 
presented  in  this  article,  tne  results  from  an  application  of  free 
vortical  system  methods  to  calculate  rotor  wake  and  the  induced 
velocities  caused  by  it  in  nap-of-the-ear th  flight  demonstrate  tnat 
the  metnods  in  this  article  for  calculating  tne  strength  and  position 
of  ground  vortices  can  be  applied  to  an  analysis  of  the  aerodynamic 
characteristics  of  rotors. 

At  many  places  in  tnis  article,  we  have  quoted  tne  experimentally 
measured  data  from  reference  L  5  J  in  order  to  confirm  the  calculation 
results  in  tnis  article.  Because  of  tnis,  we  wish  to  give  our  special 
thanks  to  Professor  H.C.  Curtiss,  Jr.  and  Dr.  Sun  Mao. 
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A  STUDY  OF  THE  ROTOR  WAKE  IN  NAP-OF-THE-EARTH 

He  Cheng  jan  and  Gao  Zheng 
( Nanjing  Aeronautical  Institute ) 

Abstract 

In  this  paper  an  investigation  of  the  ground  vortex,  the  significant 
aerodynamic  phenomenon  in  rotor  aerodynamics  on  nap-of-the  earth,  is 
carried  out.  Based  on  the  analysis  of  the  rotor  wake  near  the  ground,  a 
theoretical  method,  including  the  analytical  model,  the  formulae  and  the 
computer  program,  has  been  established  which  can  be  used  for  calculation 

of  the  ground  vortex - the  longitudind  position,  the  height  above  the 

ground  and  its  strength.  By  using  a  free  wake  model,  with  the  ground 
vortex  taken  into  account,  the  distribution  of  induced  inflow  at  the  rotor 
disc  is  also  calculated.  All  computational  results  are  compared  with  avai¬ 
lable  experimental  data  and  found  in  good  agreement. 


THE  SEPARATION  CRITERION  FOR  STEADY 
3-DIMENSIONAL  VISCOUS  FLOWS 


Xueying,  and  Liu  Mouji 


ronautics  and  Astronautics) 


SUMMARY  This  article  takes  the  topological  principles  of 
singularity  point  theory  as  a  starting  point  and  explains  the 
existance  of  half  saddle  points  on  the  cross  section  flow  lines 
perpendicular  to  the  lines  of  separation  and  the  cross  section  of  the 
object.  It  goes  on  to  explain  that  the  perception  by  such  people  as 
Lighthill  that  extreme  limit  flow  lines  must  raise  up  before  the 
separation  flow  lines  is  not  reasonable.  The  reasonable  way  of 
explaining  the  phenomena  should  be  that  the  extreme  limit  flow  lines 
are  always  coincident  with  shear  stress  lines.  It  is  only  after 
entering  points  of  singularity  that  extreme  limit  flow  lines  finally 
raise  up  and  leave  the  surface  of  the  object.  Finally,  this  article 
demonstrates  that  the  simpleminded  application  of  two-dimensional 
separation  criteria  to  three-dimensional  separation  situations  is 
inadequate.  Moreover,  it  presents  formulae  to  express 
three-dimensional  separation  criteria. 


Introduction 


Following  increases  in  the  flight  speed  of  aircraft  missiles,  in 
order  to  appropriately  respond  to  the  requirement  to  raise  the 
lift-drag  ratio  for  transonic  and  supersonic  speeds,  it  is  necessary 
to  select  for  use  thin  aircraft  wings  having  large  back  sweeps  and 
small  length  to  width  ratios.  The  fuselage  used  is  also  relatively 
long  and  thin  in  order  to  facilitate  the  lowering  of  wave  drag  as  well 
as  trim  drag.  With  thin  wings  and  large  angles  of  back  sweep  it  is 
very  easy  to  create  three  dimensional  separation.  Besides  this, 
because  of  the  fact  that  the  capabilities  of  fighters  and  missiles 
have  increased  without  ceasing,  it  is  often  necessary,  with  even 
larger  windward  angles,  to  carry  out  flight  maneuvers,  leading  to 
strong  three-dimensional  separation  flows  appearing  on  the  back  wind 
surfaces  of  the  aircraft.  Because  of  this,  research  concerning 
three-dimensional  separation  has  become  a  totally  obligatory 
requirement . 
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Research  into  three-dimensional  separation  has  only  been  done  in 
the  last  twenty  or  thirty  years.  However,  at  its  earliest,  it  can  be 
traced  back  to  1948,  when  Sears,  in  an  article  on  the  "boundary  layers 
of  drift  cylinders"  presented  the  concept  of  what  he  called  extreme 
limit  flow  lines.  He  recognized  that  adhering  closely  to  the  surface 
of  objects  there  were  flow  lines  existing  which  were  like  a  very  thin 
sleeve  around  the  object.  Beginning  from  the  1950's,  corresponding  to 
the  development  of  aircraft  missiles  at  that  time,  problems  of  three 
dimensional  separation  began  to  attract  the  serious  attention  of 
people.  E.C.  Maskell,  on  the  basis  of  extensive  experiments  with  oil 
flow,  drew  together  the  experimental  results  of  other  people  at  the 
time  and  presented  a  general  definition  of  three  dimensional 
separation  (l).  Such  men  as  R.  Legendre  (1956)  and  M.J.  Lighthill 
(1963).  basing  themselves  on  flow  line  equations  and  shear  stress  line 
equations,  applied  the  qualitative  theory  of  differential  equations 
and  presented  the  singularity  point  topologic  method  of  analysis. 

On  the  basis  of  this  point  of  view,  the  multivarious  flow  states  on 
walls  were  determined  by  nothing  more  than  a  limited  point  of 
singularity  from  flow  line  equations  or  surface  shear  stress  line 
equations  and  the  equations  mutually  interconnecting  tnem.  Tnese 
points  of  singularity  all  had  clear  physical  significance.  According 
to  the  rules  of  topology,  two  flow  states  which  have  different 
configurations,  if  they  have  the  same  number  of  types  of  points  of 
singularity  and  connecting  equations,  are  the  same  and  the  two  flow 
states  are  topologically  equivalent.  Because  of  this,  utilization  of 
singularity  point  topological  analysis  makes  it  possible  to  interpret 
the  significance  of  the  appearance  of  points  of  singularity  in  flow 
configurations,  to  analyze  flow  configurations,  and  to  predict 
possible  new  flow  configurations  [2]. 


1  9 


L  A  AAA  -Vuv  --- 


,•  V*  v  V  ’ . 


, 

*  V. 


Peake,  Tobak  (2),  Dallman  (1980)  (4)  and  K.C.  Wang,  in  recent 
years,  have  utilized  the  topological  concept  of  points  of  singularity 
to  analyze  many  types  of  flow  configurations  on  aircraft,  forming  one 
important  part  of  the  contents  of  the  qualitative  research  on 
three-dimensional  separation.  Due  to  the  introduction  of  topological 
points  of  similarity,  the  flow  configuration  analysis  based  on 
experimental  observations,  which  had  been  used  in  the  past,  developed 
in  the  direction  of  the  utilization  of  mathematical  tools.  Even  to 
the  present  day,  the  use  of  numerical  value  methods  to  solve  for 
three-dimensional  separation  flows  is  still  extremely  difficult. 
Because  of  this,  qualitative  research  is  still  a  very  improtant  path 
to  be  pursued.  Through  it,  it  is  possible  to  understand  the  structure 
of  separation  flows.  It  is  also  possible  to  make  progress  toward 
presenting  a  simplified  mathematical  model  to  serve  in  numerical  value 
calculations.  This  article  concentrates  on  the  discussion  of  the  form 
and  criteria  of  three-dimensional  separation. 


I.  Two  Different  Types  of  Models  in  Research  on  Three-Dimensional 
Separation 

In  order  to  understand  the  special  points  of  three-dimensional 
separation,  we  will  first  reconsider  for  a  moment  the  situation  in 
two-dimensional  steady  state  separation.  Two-dimensional  steady  state 
separation  is  a  type  of  singularity  point  separation.  The  point  of 
separation  is  located  at  the  intersection  of  flow  lines,  showing  the 
formation  of  half  saddle  points  (Fig.  1).  The  wall  surface  flow  lines 
flow  into  the  point  of  separation  from  two  sides.  After  they  flow 
together,  the  flow  lines  and  the  wall  surfaces  form  an  open  wall 
surface  with  a  certain  included  angle.  According  to  the  definition 
for  points  of  singularity,  and,  given  the  condition  that  the  normal 
velocity  near  separation  points  is  not  zero,  it  is  possible  to  derive 
the  separation  criteria  presented  below  13, 5J. 
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In  these,  *  runs  along  the  surface  of  the  object.  *  is  a  direction 
normal  to  th>  surface  of  the  object.  is  the  velocity  component  in 

direction  *  (Fig.  1). 
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Fig.  1  A  Diagram  of  Two-Dimensional  Separation 

Fig.  2  Maskell  Three-Dimensional  Separation  Model  1 .  Separation 
Surface  2.  Separation  Lines  3-  Object  Surface 

As  far  as  three-dimensional  separation  is  concerned,  there  are 
always  two  different  points  of  view  (models)  in  existence,  that  is, 

1.  The  Maskell  three-dimensional  model  11].  In  1955,  in  his 
famous  essay  "Three-Dimensional  Separation  Flow,"  he  first  presented 
the  idea  that  separation  lines  are  an  envelope  of  extreme  limit  flow 
lines  (Fig.  2).  The  extreme  limit  flow  lines  flow  toward  the 
separation  lines  from  the  two  sides  of  the  separation  lines. 
Moreover,  the  same  separation  flow  line  cuts  itself.  After  this,  it 
leaves  the  surface  of  the  object  and  forms  a  three-dimensional 
separation  flow  surface  (it  is  a  one  flow  surface.)  Separation 
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surfaces  come  in  two  forms,  three-dimensional  separation  gas  bubbles 
and  free  vortical  surfaces  (Fig.  3)*  The  gas  bubble-form  separation 
line  normally  connects  saddle  point-form  points  of  singularity.  At 
points  of  singularity , =  This  is 

equivalent  to  the  disappearance  of  shear  stress,  that  is,  t=t=0. 
On  separation  lines,  besides  this  type  of  independent  point  of 
singularity,  they  are  all  ordinary  points.  At  the  locations  of 

jj  dV 

ordinary  points  and  -jt  cannot  both  be  zero  at  the  same 

time,  that  is,  the  shear  stress  cannot  be  zero.  Therefore,  at 
ordinary  point  locations,  one  has  at  least  two  extreme  limit  flow 
lines  mutually  intersecting  to  form  a  peak  point.  The  separation 
lines  are  the  tracks  of  this  type  of  peak  point. 


Fig.  3  Two  Different  Types  of  Separation  Forms  (a)  Gas  Bubble  (b) 
Free  Vortical  Surface  1.  Viscous  Layer  2.  Gas  Bubble  Surface  3* 
Separation  Line  4.  Extreme  Limit  Flow  Lines  5*  Surface  of  Object  6. 
Free  Vortical  Surface  7.  Flow  Lines 


According  to  Maskell's  definition  of  separation  lines  and 
separation  surfaces,  it  is  possible  to  derive  the  fact  that  the  cross 
section  flow  lines  on  the  cross  section  plane  perpendicular  to 
separation  lines  and  the  object  surface  intersect  each  other  at  points 
on  the  flow  lines  and  cross  sections  to  form  semi-saddle  points  (Fig. 
4">.  This  is  similar  to  the  situation  with  two-dimensional  separation. 
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The  cross  section  of  extreme  limit  flow  lines  with  the  plane  xoz 
forms  a  point  of  singularity  at  point  .  When  there  is  separation, 
it  is  necessary  to  have  a  normal  velocity  to  leave  the  surface  of  th 
object.  Because  of  this,  equation  (1 )  is  satisfied  at  point  o}  . 

Flow  lines  on  the  separation  surface  cut  across  separation  lines. 

They  do  not  cut  across  the  surface  of  the  object.  Because  of  this, 
the  separation  surface  forms  an  angle  of  a  certain  size  with  the 
surface  of  the  object. 


Fig.  4  Maskell  Separation  Model  and  Cross  Section  Flow  Lines  y—  - 
Separation  Line  Direction,  -  Direction  Normal  to  the  Surface  of 

the  Object,  *-  -Positioned  on  the  Surface  of  the  Object  Perpendicular 
to  >• 

Fig.  5  Separation  Lines  on  Two-Dimensional  Straight  Aircraft  Wings 
and  Sweptback  Wings 


Obviously,  Maskell' s  concept  of  three-dimensional  separation  can 
be  drawn  out  as  an  extension  of  two-dimensional  separation.  If  we 
assume  a  two-dimensional  straight  aircraft  wing  having  on  it  a 
separation  line  parallel  to  the  forward  edge,  when  the  forward  edge  of 
the  wing  is  swept  back,  tne  two-d imensional  separation  line  then 
becomes  Maskell' s  three-dimensional  separation  line.  It  is  equivalent 
to  two-dimensional  separation  with  the  addition  of  a  component  for  the 
direction  of  development  (Fig.  5). 
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Maskell's  separation  model  finds  important  application  in 
research  on  three-dimensional  separation.  In  recent  decades  it  has 
been  extensively  quoted  by  people.  In  oil  flow  spectrum  analysis,  he 
presented  methods  for  determining  lines  of  separation  on  object 
surfaces  as  well  as  auxilliary  lines.  Moreover,  he  was  able  to 
connect  the  flow  configurations  for  object  surfaces  to  the  flow 
configurations  in  the  air. 

2.  The  Lighthill  and  Peake-Tobak  model.  After  the  1960's,  people 
became  more  and  more  interested  in  the  flow  configurations  for 
aircraft  at  large  windward  angles.  At  such  times,  leeward  surfaces 
often  show  the  appearance  of  very  complicated  separation  flows, 
including  vortical  flow  forms  and  point  vortices.  According  to 
Maskell's  model,  there  was  no  way  to  solve  for  these.  Because  of 
this,  these  types  of  point  vortices  were  often  seen  as  oil  vortices 
and  were  considered  to  be  phenomena  produced  by  oil  flows. 

Observations  of  air  flow  configurations  demonstrate  that  point 
vortices  on  the  surfaces  of  objects (also  a  type  of  separation 
structure  point)  are  start  points  for  simple  root  filament  vortices, 
and  are  often  the  start  points  for  the  vortical  cores  of 
three-dimensional  separation  vortices  (Fig.  6). 

In  the  Maskell  model,  extreme  limit  flow  lines  form  peak  points 

at  points  of  separation.  At  one  point,  there  are,  simultaneously,  two 

upward  flow  lines  crossing.  Moreover,  the  lines  of  separation  cut 

each  other.  Lighthill  [2]  demonstrated  that,  in  the  vicinity  of  lines 

of  separation,  it  is  not  possible  for  extreme  limit  flow  lines  to 

mutually  intersect  separation  lines-that  they  lift  off  the  surface  of 

objects  in  front  of  the  lines  of  separation.  Because  of  this,  such 
people  as  Lighthill  and  Peake-Tobak  did  not  use  extreme  limit  flow 

lines  to  define  and  analyze  the  flow  characteristics  in  the  vicinity 

of  three-dimensional  separation  lines.  They  made  use  of  shear  stress 

lines  to  make  their  definitions  and  analysis.  They  recognize  that 

shear  stress  lines  exist  everywhere  on  tne  surface  of  objects  and  that 

shear  stress  fields  on  the  surface  of  objects  are  continuous  vector 

quantity  fields.  The  extreme  limit  flow  line  equation  is 
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,  from  non-glide  conditions  on  the  surface  of 
p-*0  ,  «*-*o  and  L’Hopital's  law 
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From  half  ton  f  ormulas  az  ’  \  A  #, 

Therefore,  equation  (3)  can  be  written  as 


dy/dx-  t,/t. 


The  shear  stress  line  equation  (4)  and  the  extreme  flow  line  equation 
(3)  correspond  completely.  In  the  vicinity  of  lines  of  separation, 
extreme  limit  flow  line  equations  do  not  represent  actual  flow  lines. 
They  simply  represent  their  projections  (this  point  will  be  discussed 
later.)  Because  of  this,  Lighthill  only  uses  equation  (4)  to  discuss 
separation. 
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1  Lighthill 


Fig.  6  Point  Vortices  on  the  Surface  of  Objects  1 .  Point  Vortex 
Fig.  7  Lighthill  Separation  Model  2.  Shear  Stress  Line  3-  Extreme 
Limit  Flow  Line  4-  Object  Surface  25 
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Starting  from  the  point  of  view  of  a  topological  analysis  of 
points  of  singularity,  Lightnill  and  Peake-Tobak  recognize  that  ines 
of  separation  arj  basically  shear  stress  lines.  The  starting  points 
for  lines  of  separation  are  saddle  points,  and  their  end  points  are 
the  end  points  of  the  separation.  On  the  two  sides  of  the  separation 
line,  the  shear  stress  lines  take  the  separation  line  as  th°ir 
asymptote  leading  to  convergence  .  They  finally  flow  together  at  the 
end  point  of  the  separation  (Fig.  7). 

7.  The  debate  concerning  the  definition  of  separation  lines.  The 
debate  over  whether  lines  of  separation  are  asymptotes  or  envelopes  is 
still  continuing  up  to  the  present.  Because  of  the  fact  that,  in  such 
experiments  as  those  involving  oil  flows,  the  resolution  in  the 
vicinity  of  the  separation  lines  is  limited,  it  cannot  be  determined 
whether  they  are  asymptotic  or  they  cross  Mathematically,  it  is  at 
present  not  possible  to  directly  solve  the  N-S  equations  involved  in 
three-dimensional  separation  problems.  Inside  China,  Znang 
Hanxin-bas ing  himself  on  a  certain  separation  model-maintai ns  that,  as 
far  as  boundary  layer  equations  are  concerned,  separation  lines  are  an 
envelope,  but,  in  terms  of  N-S  equations,  separation  lines  are 
asymptotes.  Moreover,  he  recognizes  that,  on  a  cross  section 
perpendicular  t<">  the  separation  lines,  cross  section  flow  lines,  at 
points  of  intersection  between  separation  lines  and  cross  sections, 
should  be  saddle  points.  Therefore,  he  takes  two-dimensional  criteria 
concerning  separation  and  uses  them  on  three-dimensional  separation 
flows  [  6  ]  . 


II.  A  Discussion  of  Singularity  Point  Characteristics  on  Cross 
Sections  Perpendicular  to  Three-Dimensional  Separation  Lines 


The  models  of  people  like  Lightnill  and  the  MasK-lx  mode^  'in¬ 
different.  They  maintain  that  the  extreme  limit  flow  lines  cur.  up 
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in  t:.r  vicinity  ol  separation  *  ir.es  ana  never  interse  ct  s-p  irv.  .  . 

1  in-*.--  at  any  pc.  nr.  T.*.ey  u  not  cv-n  intersect  eacr.  utr.rr.  in  tn 
vicinity  of  tne  separation  line  pro  jeot  i  ions ,  tney  Itmow  tne 
direction  of  tne  surface  shear  stress  lines  and  f ina^.y  flow  tc-getr.er 
at  tne  end  point  of  tne  separation.  if  tney  are  perpendicular  *  .  it- 
separation  lines  and  the  plane  of  the  surface  of  the  object,  tnen,  tne 
flow  line  orientation  of  that  plane  is  not  able  to  fori  na* f  saddle 
points  (Fig.  i : .  Because  of  this,  in  topological  metnoas  of  analysis 
for  separation  points,  it  is  necessary  to  carry  out  adjustments  of  tne 
relationships  of  extreme  limit  flow  lines  and  shear  str-ss  lines  in 
the  vicinity  of  separation  lines  in  order  to  make  tr.ex  reasonable . 

B  8  Lifhthill  B  9 

Fig.  3  In  the  Model  of  Such  People  as  Lignthiil,  Cross  Section  Flow- 
Lines  Are  Not  Able  to  Form  Half  Saddle  Points 

Fig.  9  Delta  Wing  Cross  Section  Flow  Line  Flow  Configurations  AlCacen 
to  Body 

1.  Points  of  singularity  perpendicular  to  tne  plane  of  separation 
flow  1  i  r.e? .  X  a. -Kell  is  different  from  tne  topological  analysis  mode* 
of  s<-par  it  1  or.,  and  the  two  of  tr.em  cannot  be  confused.  In  a 
topological  analysis  of  points  of  singularity,  one  s'r.oul  a  make  t:.- 
prints  of  separat ion  on  flow  line  cross  sect i 1 ons  points  f 
s  :  r.  r*  lari  ty-half  saddle  point?. 

A  r.  /•  ♦  t  K.  t.  r<  a  ->  r  v  f  tr  t  r  n r  '  :  •  a  .  a  r,  '  v  s  ’  -  ■  f  r  a  i  n  *  f 


From  to  r-x-iaj.  *e.-  «?.  v-.n  above ,  or.tr  jar.  s-e  tt.at  ,  i  r.  a  t : p 
analysis  of  points  of  similarity,  i  ntersect  1  or.  points  Or  tween 
separation  lines  and  miss  sect  ions  are  half  said.-  p,.r.ts.  * 
being  the  case,  now  one  understands  tr.c  struct  u  e  of  separa*  , 
surf  a  'es  and  tn~  format  i  or.  of  naif  sadd. e  pc  .  nt  s  i  s  a  quest  .  or, 
of  discussion. 

2.  T  r .  forma*  i  or,  of  se  para*  ion  surface;--  in  toj  •.  .of:  a.  ar, 
of  po  i  r.  *  s  of  singularity.  In  the  VasK-i.!  mode.,  separat .  or.  p. 
forme  j  from  ex*reme  limit  flow  lines  s-parat  i  ng  from  tr.-  c.rfa 
objects.  In  tne  Ligntni  11  model,  separation  surfaces  are  form 
the  turning  up  of  extreme  iu.it  f.  ow  lines.  They  sn.  u.  i  o.jm>- 
or  igination  point  of  upper  flow  po.r.ts  of  s.nguiar .  ty  ai.u  sej  I 


F:  rf.  1  1  ,  a. .  pjintj  of  singular .  ty  on  trie  surf-i.v 
'"i  r  He?  p  o  *  r*  t  s  ot  Sirif'^ar.  ty  in  tu  r  l  x  j*  .  *  i » * 
one  see;  trie  uppcur-iii.'e  of  pjiiitd  vi  singular. ty 
.v  p  e  r  pr  na  *  •  a.  nr  p.-tins  !"!.  . .  s  io  o!<v,-^  .y 

*  taut:  or.  we  dismissed  before  in  w.n.en  oruy  ;r  osc 
i'-i*  :  ons  r.'jd  points  of  s  i  ngul  t»r  i  t  y  ana  trie  flow 
:.e  otr.er  t  w  per  per.  a  i  tuiar  pi  a  r.  e .  ■  aid  n  o  t  .  For 
stationary  point,  in  term.-  of  t :. e  oejeet  surface 
point.  Wail  purfnt*  sr.ear  stress  i  mes  are  a. so 
**•  u’uhI  iy  ortr.ogonai  pi  ur.e?  wr..o:. 
are  perpend,  tu.ar  t  _  tr.e  .surf  a  .•  e  _ f 
*  r.e  ex.s'ar.o-  of  :  oil  f  sadd.-  p  .r.*.  . 
t  -p  1  '-a:  'a.  ar.alvs  .  f  p  .  .  r. *s  .  f 

f.oW  *  1  fir  s 


is  f^ried  from 


W  ? .  - 


of  • 


As  is  shown  ir 
an  orthogonal  :urv- 
of  th"  object .  C 


3  u  r  f  a 


o  v ■ 


?:*.  1 ^ ,  on  the  surface 
coordinate  system.  fcon 
is  norma.  * 

t-0. 


•.  ot  ject ,  we  sele  ?t 
.  :  wiel  to  the  surface 

surface  ul  the  object.  on  trie 
Tne  extreme  flow  line  equation  or 


surface  snear  stress  line  equation  is 


Ml/'.-Mi/’., 

V  t> 


A,  .  j  Aj  art  lam-  coefficients. 


O.-.. 


-v-.y  par  a.  .el 

4o  tr.r  component 

o  t  ^  .  • 

n.  * 

a.-.-  x"-  a  *■  *  - 

.  s  4 1-  d  :  s  t  an  ;-e 

between  extreme 

1  im .  t  f  * ...  w 

a  *  lie s  .  A 

.  *  r."  e  .  gr  t 

:f  the  flow  tube. 

Tne  amount  of 

mass  flow.n 

,g  tf.TOag 

w  *  u  t  - 


"l  *■  p/irui 


p  .  :  s  dens  :  *v .  *  :s  average  v~.  .  • :  *  ,v  .  l-rfa  > 


in 


Because  of  tnis 

«-A*nx./2 -const 

*-c(  »?r)  "•  (,0) 

when  approaching  separation  lines,  ■  .  Because  of 

thir,  in  the  vicinity  of  three-dimensional  separation  lines,  the  flow 
lines  inside  flow  tubes  must  turn  up.  This  includes  extreme  limit 
f  1  o  w  1  i  n  e  s  . 

T.ne  proof  described  above  explains  that  the  raising  of  flow  tubes 
in  tne  vicinity  of  separation  lines  is  correct.  However,  it 
recognizes  that  the  separating  of  extreme  lir*it  flow  lines  from  the 
surfac--  of  objects  is  not  certain.  If  it  follows  the  movement  upward, 
tr.-r.  its  lower  surface  must  also  have  flow  lines.  This  is 
‘ont  rail  story  to  tne  definition  of  extreme  limit  flow  lines.  Because 
of  *  n  i  s ,  extreme  limit  flow  lines  will  not  leave  the  surface  of  the 
ob-ect.  i*  -  only  at  tne  location  of  points  of  singularity  that 
ex*r-m-  limit  flow  lines  flow  togetner.  After  tnat,  they  separate 
from  tr.~  surface  of  t.ne  object  at  a  given  angle.  Fig.  14  presents  the 
p-.  *  *  -  ms  for  extreme  limit  flow  line  separation  from  the  surface  of 
-  r  j  -  •  *  s  in  t.ne  Yask-11,  Ligr.thill  models  as  wel.  as  for  the  model 


?  .g.  1  4  :  .  .  js*  ra*  .  -.-ns  of  the  Alignment  of  Extr-m* 

~,n  r--  i  .  ff-r-r.*  Typ-s  ■  >f  Yodels  1  .  Ma.-sitei  i  hod  - . 

7 .  Tn-  Kod~l  In  Tn.s  Article 


Li  mi  *  Fiow  n-s  in 
y  .  s. o d (  x 
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4.  A  discussion  concerning  upward  flow  line  singularity  points 
perpendicular  to  separation  line  cross  sections.  Due  to  the  fact  that 
separation  planes  are  formed  from  flow  lines  which  come  out  from 
points  of  singularity  which  are  nodal  points  that  appear  on  cross 
sections  perpendicular  to  separation  lines,  when  they  intersect  cross 
sections  (sections  above  the  surface  of  objects)  perpendicular  to 
separation  lines,  it  forms  a  break  line  (the  os  line  in  Pig.  12(b)). 
Due  to  the  fact  that  extreme  limit  flow  lines,  on  lines  of 
separation,  do  not  turn  up  again,  on  cross  sections  perpendicular  to 
separation  lines,  they  form  respectively  straight  lines  from  the  two 
sides  pointing  toward  intersection  points,  for  example,  and  So 

in  Pig.  12(b),  which  give  rise  to  the  appearance  of  half  saddle 
points . 


III.  Conditions  Which  Should  Be  Satisfied  in  the  Vicinity  of 
Separation  Lines 

We  have  already  explained  that  separation  lines  are  both  surface 
shear  lines  and  extreme  limit  flow  lines.  If  we  take  the  separation 
line  as  the  y  axis  with  the  *  axis  perpendicular  to  the  separation 
line  (separation  line  coordinate  system),  the  corresponding  shear  - 

force  components  are  .  On  the  separation  lines,  t,#  o  ,  t,  ,=  o  . 

Alonq  the  x  direction,  the  shear  force  component  *,  changes 
sign  in  the  vicinity  of  separation  lines.  Because  of  this,  ax,/ax<  0 , 

It  follows  from  this  that 


T«“0,  »1./dx<0 


(11) 


This  is  a  necessary  condition  for  three-dimensional  separation. 
However,  the  condition  described  above  •does  not  consider  the  influence 
of  changes  in  shear  stresses  along  the  direction  of  lines  of 
separat i on . 
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I /.V.v.  ^V.vlvJV- 


Consider  for  a  moment  the  situation  of  normal  velocity  u'  in  the 
vicinity  of  separation  lines.  In  order  to  do  this,  let  us  take  a 
certain  point  on  >\  separation  line.  Let  us  also  take  a  point  •#» 

which  is  at  a  distance  from  the  separation  line.  For  convenience 

sake  let  us  choose  point  p ,  as  the  origin  point.  Then,  the  normal 
velocity  at  point  R  can  be  represented  as 


»--»'+(" ), 


(12) 


From  the  continuity  equation 


~dX~  +  ~iy~  +  0  th3t  1S 


»Z  \  dX  +  iy  J 


(13) 


On  wall  surfaces  we  have 


0  , 


Because  of  this 


(*u>/dz')r=  -t(*u/dx),+  (dv/iy),]  =  0 


Due  to  this  (#*u>/i7*)P=  —  Z(*tu/dzdx'),+(dlv/dydz),') 


Wgm. 


1 

2 


-C(dt*/azdx')r+(  alo  I  »yiz  )  3  a* 


X 


Therefore 


(H) 


It  is  required  that,  at  the  time  of  separation, 
of  this,  it  is  required  that 


v>r>  0  , 


Because 


(dtu/dZdx)f+(*lv/dyd?)f<  0 


That  is,  for  each  point  on  the  separation  line,  one  should  have 


(at./dx)+(it ,/dy)<  0 


(16) 


One  can  see  from  this  that,  in  the  vicinity  of  three-dimensional 
separation  lines,  the  conditions  which  should  he  satisfied  are 


0  , 


Un¬ 

ix 


(17) 


Reference  [6],  in  solving  for  saddle  points  on  cross  sections,  makes 
use  of  the  condition  tu,>o  .  However,  it  does  not  take 

equation  (16)  and  substitute  into  it  the  conditions  which  should  be 
satisfied  for  three-dimensional  separation.  Below,  we  discuss  the 
necessity  of  satisfying  the  condition  (*T,/*x)  +  (#T„/*y)<  o 

As  far  as  closed  type  separation  is  concerned,  that  is,  the 
situation  in  which  the  point  of  origination  is  a  saddle  point  and  the 
end  point  is  the  separation  nodal  point,  we  still  select  a  separation 
line  coordinate  system.  Moreover,  we  assume  that  A,  Bt 
respectively  are  the  initial  and  final  points  of  singularity  (Fig. 

15).  If  there  is  a  fir  w  line  through  a,  B  ,  and  the  direction  of 
that  flow  line  is  from  to  B,  ,  then,  on  this  flow  line,  changes  in 
the  rules  for  t„ ix/ty  for  various  points  should  be  as  shown 

in  Fig.  15-  The  line  of  separation  is  just  this  type  of  flow  line. 

As  far  as  *T»/*y>  0  j_s  concerned,  this  section  must  guar^ni-pp  (^T./*y)  +  (*T,/*.vX  0 

It  is  then  required  that  #*./#*<—  ixjty . 

That  is,  it  is  not  simply  required  that  it./ix<  0. 


Pig.  15  Distribution  of  T>  Along  Separation  Line 


Conclusion 

As  far  as  research  on  the  criteria  for  three-dimensional 
separation  goes,  there  is  still  a  good  deal  of  work  that  needs  to  be 
done.  For  example,  the  origination  point  for  open  form  separation  and 
its  criteria.  Even  if  one  is  dealing  with  closed  form  separation,  the 
sufficient  conditions  for  it  have  not  yet  been  obtained.  Following 
along  with  the  development  of  calculations  and  experimental 
techniques,  it  is  necessary  to  reach  a  deeper  understanding  among 
people  of  the  phenomena  produced  by  three-dimensional  separation. 

The  various  points  of  view  in  this  article  were  discussed  with 
such  comrades  as  Zhang  Hanxin  and  Wu  Jiezhi,  from  which  we  obtained 
not  a  little  instruction  and  profit.  For  this  we  express  our  thanks. 
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STUDY  OF  THE  SEPARATION  CRITERION  FOR  STEADY 
S-DIMENSIONAL  VISCOUS  FLOWS 
Lii  Zhiyong,  Deng  Xueying  and  Liu  Mouji 

( Beijing  Institute  of  Aeronautics  and  Astronautics ) 

Abstract 

In  the  case  of  steady,  2-dimensional  separation,  separation  point  is  of 
a  critical  type.  The  limiting  streamlines  on  the  wall  from  both  sides  of  the 
separation  point  run  into  the  point  and  then  depart  away  from  the  surface 
of  the  body  with  a  certain  angle.  According  to  the  definition  of  a  critical 
point  and  the  condition  that  the  normal  velocity  closed  to  the  separation 
point  is  non-zero,  the  separation  criterion  for  steady,  2-dimensional  flows 
can  be  deduced  as  followsi  The  friction  stress  at  the  separation  point 
vanishes  and  There  is  a  reversed  flow  after  the  separation.  In  the  case  of 
steady  3-dimensional  separation,  there  is  a  separation  line  which  starts  at 
a  saddle  point  and  ends  at  a  separated  nodal  point,  as  well  as  a  separation 
sheet  which  consists  of  streamlines  originating  from  a  critical  point.  This 
critical  point  can  be  viewed  as  a  saddle  point  topographically  on  the  wall, 
or  as  an  attached  nodal  point  in  the  section  containing  the  separation 
line  and  perpendicular  to  the  wall.  Based  on  the  theory  of  singularities 
anfl  topological  rule,  there  may  exist  a  half  saddle  point  in  the  section 
cut  normally  to  the  wall  and  separation  line.  As  a  result,  it  can  be  argued 
that  the  argument  about  the  limiting  streamlines  must  rise  before  approa¬ 
ching  separation  line,  due  to  Lighthill,  is  incorrect.  This  paper  shows 
that  the  limiting  streamlines  along  the  wall  always  coincide  with  the 
corresponding  friction  lines.  Limiting  streamlines  are  attached  to  the 
wall  until  the  separated  nodal  point  or  saddle  point  where  they  lift  off 
the  wall.  Finally,  the  paper  also  presents  the  separation  criterion  for  steady 
3-dimensional  flows  which  is  different  from  that  for  2~dimensiona1  flows. 
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DETERMINATION  OE  OPTIMAL  POSITION  0 
FOR  FLEXIBLE  FLIGHT  VEHICL 


Yuan  Jianping  and  Chen  Shi 


SUMMARY  This  article  does  research  on  the  determination  of 
optimal  positions  for  control  mechanisms  in  flexible  aircraft.  In 
control  systems  which  have  rate  of  speed  feedback,  the  control 
function  has  been  seen  as  an  equivalent  damping,  and  w^  have  derived 
the  criteria  for  optimal  positions  of  the  rudaer  face  as  determined  oy 
aerodynamic  elasticity  effects.  We  nave  also  presented  formulas  for 
adjusting  the  relationships  between  the  position  of  rudder  surfaces 
and  the  position  of  gyroscopes.  Practical  examples  are  appended  for 
empirical  demonstration. 


I.  Introduction 


Control  mechanisms  refer  to  control  surfaces  or  oscillating 
engines.  They  can  be  abstractly  conceived  of  as  components  producing 
concentrated  forces  or  moments  of  force.  The  effects  of  control 
mechanisms  are  an  important  source  of  disturbance  for  the  elastic 
vibration  produced  in  aircraft.  It  can  give  rise  to  curvature  or 
twisting  vibration  of  aircraft  structures.  Therefore,  it  changes  the 
size  and  distribution  of  aerodynamic  forces  in  aircraft.  It  alters 
the  efficiency  of  control  surfaces  and  influences  the  control  systems 
of  aircraft.  Differing  installations  of  control  mechanisms  will 
produce  differing  results  in  the  excitation  of  vibrations  as  far  as 
elastic  vibration  movements  are  concerned.  Because  of  this,  the 
selection  of  optimum  positions  for  control  mechanisms  has  practical 
significance . 

The  purposes  of  making  the  positions  of  control  mechanisms 
optimal  are  to  be  able,  for  general  requirements,  to  reiiablly 
control  aircraft  movements  and  restrain  elastic  vibration.  Such 
people  as  G-evarter  [1]  and  Chen  Shilu  [2]  did  research  on  the 
influence  of  control  mechanism  positions  and  sensor  positions  on  the 
stability  of  flexible  aircraft.  They  made  use  of  transmission 
functions  and  zero-polar  methods  respectively  and  discussed  the 
conditions  which  should  be  satisfied  by  the  locations  of  control 
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mechanisms  and  sensors.  As  a  type  of  drive  cor.t  r*,*  meti.ou,  s  4 
people  as  Arbel  [3]  and  Schultz  14]  presented  an  appropria:-  • 
of  location  for  control  mechanisms  in  order  to  cause  fc-dta.v.  sy.;  ir-mo 
to  give  rise  to  increases  in  flexible  structure  damping  effects  o,,  xs 
to  facilitate  the  production  by  limited  control  of  maximum  r-strai  n.  r.g 
effects  on  elastic  vibrations.  However,  all  of  th--m  wer--  d-rp-na-.n'  on 
numerical  value  calculations.  They  were  not  able  to  offer  ana ly i . 
relationships.  This  article,  basing  itself  on  tne  Dasic  concept  of 
the  laws  of  aerodynamic  energy,  on  tne  foundation  of  tne  word 
mentioned  above,  discusses  a  step  furtner  tne  question  of  tne 
optimization  of  the  location  of  elastic  aircraft  control  me cnani sms . 

In  systems  having  velocity  feedback,  we  derived  two  practical  ini-x 
functions  and  a  simple  relationship  showing  tne  adjustment  pos  1 1 1  on:  r»g 
of  rudder  surfaces  and  gyroscopes.  Moreover,  it  carries  out 
empirical  demonstrations  using  practical  examples. 


II.  Optimum  Position  of  Rudder  Surfaces 

This  article  considers  the  optimum  position  for  rudder  surfaces 
in  control  systems  with  velocity  feedback.  The  equation  describing 
the  elastic  vibration  of  aircraft,  the  measurement  equation,  and  the 
control  law  are 


q+D,q  +  K,q’*BlU  +  0 

Zt-H,q 

U--GZX 


In  these,  q  is  the  generalized  coordinate  for  elastic  vibration,  n 
is  a  dimensional  vector  quantity.  D,  is  an  "  x  *i  generalized, 
natural  damping  matrix.  is  an  *  x  * |  generaliz  'd  rigidity  matrix.  Q 

is  an  uncontrolled,  generalized  force.  ni  is  a  dimension  vector.  Z, 
is  the  output  vector  for  the  T  dimensional  sensor.  V  is 
the  *"•  dimensional  control  vector.  G  is  an  ^  x  r .  feedback  matrix. 


40 


mV *  ■  ’/I  ’ 


,s 


t.)  - 

t.)  - 


expr-s.-e.;  th-  val  u-  of  a  vi  brat :  or.  mode  derivative  of 
order  »  -i*.  t:.e  loca*.  ior.  4y  of  the  number  j  'sensor.  /r((  expresses  th 

3’r.  ■*'<  force  produ  :ed  in  order  »  vibration  when  tn-  /  *  t  n  control 
surface  -nas  a  unit  deflection.  Let  us  assume  tnat  r).  is  the 
coordinate  along  toe  vertical  axis  of  tne  ,j  th  control  surface 
pre?sure  ■* e n “ e ”  and  tha*  Pn  and  are  respectively  the 

distances  in  tne  same  direct . on  from  %  to  tne  forward  and  rear  end 


P  o 


c ont r ol  surface , 


F"--%rSS 

M  p,*<h.  .  fn(»)  „ 

erJ dx  u,  y*y 


M  an:  m,  are  respec*  i  v  ly  the  mass  of  the  aircraft,  and  the 
generalized  mass.  is  the  aerodynamic  derivative  for  tic  J-  of  the 

control  surface.  When  handling  the  dispersion  integrations  described 
above,  if  the  rn*  io  1  <-♦  w^er.  the  effective  rudder  surface  area  and  the 
characteristic  surface  area  of  the  aircraft  is  very  small,  then  it 
will  apr-roxir.ate  the  form 


c:  if.  V  t  c. 

c f  fr*.  it;:  ■  k  ir.tr  d  u 


•  •»  -  *  •  j  ■ 

r. 4  .  ? 


For  ar.y  arr  i  *  nary  v^.v. 


t.  »».  .  to- 


BfiH, - xr{BiGHl  +  H\G'B'x)+  (B,GH  t- H\(P  B\) 


Ir.  tv--  equ--. 4  i  or.  abov-,  tv.e  fir?4  ■ '  :  r  r.  c-f  4  •.•  rig:.4  r.  •. i  .  i  i-  i  ? 
symmetrical .  It  i?  the  f  ~ed i  ex  cutI’  f  i  v  I -r.i  to:. ion  enust?  t  y  ?  *  • 
to  have  an  energy  exchange-  to.tr.  ti.r  er.v  i  ronucr.t .  it  van  t  ?  utter.?  ..v 
as  the  dissip&t.on  Section.  i  I .  e  Second  iect  .  JI.  .  is  IT.  .".a  .  * 

causes  the  system  not  to  have  an  er.tr^4  exchange  w  i  t the  -  •.  v  i  r  on::.- r. 4  . 
It  can  be  understood  as  the  inertial  section  ..or  gyr  ••?  _  'tie-4.  . 

Fecause  of  the  fact  that  for  any  arbitrary  value  of  i  g,  .  4  :.e  quadrat  i  • 
fern  l3tGH HYFBl)*}  is  always  zero,  in  terms  of  its 
physical  meaning,  the  sum  of  tr.e  work  done  by  f •.•roe  (BtGH H\GTBTl)q 
in  any  infinitely  small  real  d  i  splaoemer. 4  dq  ,  is  ztr  _•.  j.-- 4  us  say 

that 


Then,  we  can  recognize  that  D,q  is  the  equivalent  damping  force 

controling  elastic  vibration  in  the  system.  In  order  to  have  the 
greatest  res  +  r-  ir.ing  effect  on  elastic  vibration  in  the  system,  we 
suggest  th-  .  *he  index  of  selection  for  the  position  of  the  ruuder 
surface  n  c-:  such  as  to  make  /  -  j  '  [«r(  t  )D.q(  t  ))<* 

J  If 


(6 


/ 
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rrrt 


Ir.  order  to  make  descr  i  pt  i  or.  ccnver.i  t.nt  ,  .<  !  u.-:  orji:  •  *:.*  '■ 

-yo  +  -r  ha?  one  velocity  gyroscope  a:  .  pur  of  rdJrr  curfa  •-  .  * 

will  only  consider  third  order  v  :  :  .  •  .  or.  form.- .  rco-.u.-r  cf  t  r. i  G 

is  a  scalar  quantity.  Wr  use  d„  •  c  represent  tnt  elem-nt  /  ; :  D,  ■ 
From  formula  [2)  to  forn.,.1-  4  ,  w~  .nav- 


•  dn  d |,  d|S ' 

D,  s  dtl  dtt  d ,, 

'  ^ii  dti  d ,, 

r  *,(  n  )<K(  4  )  V  *»  )4>j(  4 )  +  4>,(  t  )4>i(  4 1_ 

2~  _m,  ",  m« 

_  4>,(  t  )*;a  )  ,  *:(  M*;r  4  )_  -  »..(  i  )4>;(  4  )_ 

XG  ■  m,  ~“+  m,  2  «, 

♦,(_»>  )4>;r  4 1  ,  »,( *1  )4>;(  4  ),  *!)♦;(  4)  ,  4>»(  ^  )4>if  4_L 

w,  ros  "i  m» 

»,OQ»;q)  +  ».d)»K4). 

m,  m, 

iiQ  )4a(  4  )  +  ».(  n  )4»;(  4  ) 


♦,(  ^  )»;c  4 ) 


According  to  the  nature  o'  D,  .  it  is  possible  to  draw  a  few 
useful  conclusions: 

1  .  For  any  arb  •  •  r-- ".v  values  of  ’4.  n,  ,  we  have  »l=  0  • 
Pecause  of  this,  Ot  -.-t  least  has  a  Zero  character  is  t  i  e  roe  *  . 


r\  ull  un 

2.  The  fundamental  for"\  of  u>  is  ={dudtt  —  </,!)<  o  • 

If  0  ,  ,  then,  matrix  D,^  11 

useful.  It,  has  one  each  positive,  zero,  and  negative 

characteristic  roots.  When  considered  as  a  damping  fore-,  D,q  ,  in 
certain  circumstances,  creates  a  damping  effect,  and,  in  otner 
circumstances,  creates  an  accelerating  effect.  The  overall  effect  is 
determined  by  the  ratio  relationships  between  the  three  roots. 
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eq.«-  '  .  *<  -  '  '  g-  coordinate  for  undisturbed 

•.  i, ^  is  the  kinetic  energy  of  the 

*/,  .  Oi  *  •  •  g-r.*rf  1 :  zed  force.  When  we  consider  the 

. ■  o ^  *<»*•  "Ajtr'-e*' < '-o")  v  (*<  +  *<)  pn4  (*»+*<) 

*<  and  i  r.  formula  (9-.  With  just  a  little  rearranging 
j  ■  c  r.  b  rr  e 


(o.jXy  *t  0,y-*y  "t  —  i  *  1  ,  2*  ‘"i  ® 


o.y -a,..  ’  >< 


r. on- *  inear ity  terms  for  x,  and  it' 


Cn - ^(Y.y  +  Yy(),  - ^-(Y.y-Yyy). 


r.  equations  contained  in  formula  \10)  and  write  them  in 
form  as 


AX  -  -  BX  -  CX  -  E X  -  PX  +  Y 


A-(a„)  is  the  generalized  mass  matrix,  symmetrical,  fixed 


is-  the  generalized  damping  matrix,  symmetrical 


is  a  generalized  conservation  of  force  coefficient 


svmm-  r i cal 


£■“(*//) 


racial  correction  force  lu&irix,  asymmetr  i  cal 


p«(pi>.)  ’ the  gyroscopic  force  matrix,  asymmetrical 

Introducing  the  Rayleigh  dissipation  function 


In  this  case,  the  drag  force  in  formula  (11 )  can  be  written  as  aF/aX, 
TJ sine  an  orthogonal  transformation,  we  take  p  ••  nd  change  it 
into  a  square  sum  form 


I f  6t>o,  ’  "then,  force  6*-**  causes  the  movement  to  reduce 

speed.  On  the  other  hand,  if  &«<  0,  ,  then  b,xt  causes  the 

movement  to  increase  speed.  Therefore,  the  force  derived  from  F 
may  be  of  a  dissipating  form.  It  is  also  possible  for  it  to  be  of  an 
accelerating  form.  If  Xj6*>0,  ?  then,  one  can  say  that,  in  the 
system  as  a  whole,  the  dissipating  forces  are  in  a  dominant  position 
over  the  accelerating  forces. 

On  the  basis  of  the  Lyaounov  theorem,  it  is  possible  to  know 
that,  if  acceleration  forces  are  dominant  over  dissipating  forces, 
then,  the  system  described  in  formula  (10)  is  unstable.  If  the 
determinant  form  Y„|<o,  ,  then,  the  aerodynamic  system  in 

formula  (10)  is  also  always  unstable.  These  two  points  respectively 
provide  basic  principles  for  our  selection  of  velocity  feedback  and 
displacement  feedback  in  control  systems. 


I'.’H 


Because  of  the  fact  that  orthogonal  transformations  maintain 

matrices  invariable,  that  is,  2  &»  =  ,  it  follows  that 

k  i 

there  is  no  need  to  take  B  .na  change  it  into  a  diagonal  matrix.  We 
define  a  new  function  index 


/-tr.(fl)  (12) 

In  this,  tr.  expresses  the  matrix  trace  (trace).  .  pn  order  to  cause 
the  system  to  have  maximum  drag  effects  on  disturbing  movements,  the 
selection  of  system  parameters  should  cause  /  to  become  extremely 
large.  When  only  considering  gyroscope  position  and  rudder  surface 
position,  it  is  possible  to  use  D,  as  defined  by  formula  (5)  to 
replace  B,  ' n  formula  (12)-  Moreover,  in  the  simple  situation 
represented  by  formula  (5)',  it  is  possible  to  write  the  obvious  form: 


1 

m, 

_1_ 

m, 

1 

ml  , 


4>/  n  ) 
4>t(  ’l ) 

*.(»!) 


(12)’ 


III.  Actual  Examples  of  Research 


This  section  discusses  elevator  position  on  a  certain  type  of 
aircraft.  It  is  necessary  to  stress  tne  fact  that  changes  in  ruaier 
position  have  very  great  effects  on  the  rigid  control  capaDilitie  -  ol 
aircraft.  Because  of  this,  the  selection  range  for  the  paramelei  ^  ; 
extremely  limited.  Moreover,  after  selecting  it,  it  is  necessary  to 
correct  the  changes  in  rigid  aerodynamic  parameters  in  order  to 
guarantee  their  being  within  permissible  ranges. 

We  already  know  th.at  the  relevant  data  for  this  aircraft  are: 
Gyroscope  position  5  '=5. 45m;  ru.1d--r  s.-rf-ce  position  tj -6.54m. 
Vibration  form  +,(h)-0.341,  4,(  )- -0.032,  *,(  *1  )- -0.282. 

Vibration  form  derivatives  ♦!(  t  )  —  — 0.36,  4>;(U«0.07,  )-0.25. 

Generalized  mass  rn.-J4.32,  rn.-6.98,  rn.-8.48. 
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1  .  If  we  giv-'  ^  a  somewhat  different  value,  and  ti:e  other  data 
use  current  values  for  the  aircraft,  under  conditions  satisfying 
equation  ( 1  )  ,  we  calc-;'  e  the  indices  for  formula  (6).  The  result 
are  shown  in  curve  in  Fig.  1  .  It  is  possible  to  -  t-h-  whm  v' 
is  in  the  range  6.4-6.85,  all  are  possible.  When  V:=6.o  7  reaches  ar. 

extremely  large  value.  The  results  of  calculations  demonstrate  that, 
when  the  rudder  surfaces  tic  close  to  head  section  or  tail  section, 
the  system  diffuses.  This  explains  the  control  system's  continuous 
introduction  of  energy  to  the  elastic  movements. 


Fig.  1  Curves  for  Values  of  1  for  Changes  in  1  Calculated  According 
to  Different  Indices 

2.  On  the  basis  of  lormula  (7),  when  we  already  know  the 
generalized  mass  and  gyroscope  position,  the  optimum  position  for  the 
rudder  surface  is  at  the  location  where  the  vibration  form  should 
satisfy 

4h(  n  )/<M  )  =  —  10.55 

We  can  find  from  the  experimentally  determined  curves  forCMj 
that  there  are  three  values  of  **  *  ei.cn  satisfy  the  formula  above: 
1.08;  4  16;  6. 48.  ii  =  1  .Oo  is  not  permitted  by  rigid  body  control 
capabil iter .  n  =  4.16  does  not  satisfy  formula  (8).  We  eliminate 
these  two  roots.  The  results  for  ’J  =  6.48  are  in  the  ranges  obtained 
in  Fig.  1  .  Moreover,  there  is  only  a  difference  of  0.06  from  current 
values.  This  method  is  simpler  and  more  effective. 


3.  We  make  calculations  from  the  index  ( 1  2 )  formula.  niVen  a 
somewhat  different  value  for  ,  and  using  the  valuer  already  given 
for  m<,  (t)  t  the  results  are  .s  seen  in  Curve  ,£«  in  Fig.  1  . 

A*  the  place  where  n  =  6 . 4-i ,  J  reaches  its  maximum  value.  This  is 
the  same  as  the  result  from  formula  (7).  Besides  this,  at  the  place 
where  *1  =  6.54,  one  sees  the  appearance  of  a  secondary  maximum  value 
This  agrees  with  the  current  values  for  this  aircraft  very  well.  In 
calculations  other  than  these,  there  is  no  need  to  integrate  the 
movement  equations.  To  say  it  another  way,  the  index  (12)  formula 
only  depends  on  aircraft  structural  characteristics  (  mass 
distribution,  vibration  forms,  and  vibration  form  derivatives.) 

Attention  should  be  paid  to  the  existence  of  a  difference  betwee 
the  results  obtained  on  the  basis  of  the  principle  of  formula  (6)  and 
those  obtained  from  formula  (12).  The  reason  for  this  is  that  the 
calculations  from  the  index  (6)  formula  include  a  numerical  value 
integration  of  equation  set  (l).  In  it,  there  is  necessarily  a 
calculation  error.  However,  it  is  the  first  type  of  method,  as 
compared  to  the  other  two  types  of  methods,  which  requires  the  use  of 
broad  ranges  of  values. 

The  calculations  above  show  that,  starting  from  the  restraining 
of  elastic  vibrations,  the  optimum  rudder  surface  position  should  be 
at  a  distance  6.48  m  from  the  head  section.  This  position  is  moved 
forward  0.06  m  from  the  design  position  in  the  original  rigid  body 
model.  This  will  cause  the  control  moment  of  force  parameter  and  the 
drag  moment  of  force  parameter  to  go  down  respectively  by  4£>  and  1.5?° 
Other  dynamic  parameters  show  no  obvious  changes. 

IV.  Conclusion 

This  article  makes  use  of  the  concept  of  aerodynamic  energy.  It 
takes  the  effects  of  velocity  feedback  and  sees  them  as  equivalent  to 
drag.  We  have  derived  the  principles  for  optimum  rudder  surface 
locations  with  a.  view  to  restraining  elastic  vibrations.  Index  (6) 
formula  and  index  (12)  formula  are  appropriate  for  general  use  in 
velocity  feedback  control  systems.  In  the  simple  situations  which 
this  article  researches,  formulas  (7)  and  (8)  can  give  simple  and 
convenient  results.  Examples  of  practical  research  prove  the 


conclusions  above. 


methods  . 
posi t i one 


We 

Moreover,  we 
for  aircraft 


obtained  matching 
demonstrated  that 
being  studied  are 


results  from  all  three 
the  current  rudder  surface 
secondary  optimums. 
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DETERMINATION  OF  OPTIMAL  POSITION  OF  ACTUATORS 
FOR  FLEXIBLE  FLIGHT  VEHICLES 
Yuan  Jianping  and  Chen  Shilu 

(. Northwestern  Pol ytechnical  U  niversity ) 

Abstract 

Since  the  action  of  actuators  on  flexible  flight  vehicles  is  an  impor¬ 
tant  exciting  source  to  elastic  vibration,  it  is  of  practical  significance 
to  choose  appropriate  positions  for  actuators  as  well  as  for  sensors.  W e 
shall  in  this  paper  study  optimal  positions  of  rate  gyros  and  elevators 
that  are  modelled  as  lumped  operating  elements.' 


Through  the  feedback  control,  a  term  (/?,CA/,  )q  is  introduced  to  the 
system  to  count  for  the  damping  force,  where  Bi  is  a  control  matrix,  rely¬ 
ing  on  the  actuator  locations,  and  W,  is  a  measurment  matrix  depending 
on  the  gyro  locations.  An  effective  damping  matrix  D,  is  defined  in  (  5  ), 

and  then  Rayleigh's  dissipation  function  can  be  obtained  as  = 

It 

In  order  to  make  the  system  having  most  damping  effect  to  vibration, 
the  actuators  and  sensors  are  so  positioned  that  the  effective  damping 
force  (df/dq)  could  be  maximized. 

Because  of  its  real  symmetry,  D,  can  always  be  diagonalized.  In  accor¬ 
dance,  suppose  that  F  = — 2  &»<?*.  //  0  ,  the  elastic  motion  will 

i  k 

be  decelerated  by  force  if  0  ,  the  motion  will  be  decelerated  by 

force  if  b*<!  0 ,  the  motion  will  be  accelerated.  That  is,  the  forces 

derived  from  F  may  be  either  dissipating  or  accelerating.  When  6>~>  0  . 

k 

it  is  thought  that  the  dissipating  force  is  leading  to  the  accelerating  one. 

Since  the  trace  of  a  matrix  will  not  be  changed  by  a  normal  transfor¬ 
mation,  we  have  2^*  =  tr.  D,.  Then  a  new  performance  index  is  defined 
h 

by  (12).  The  problem  under  consideration  is  to  determine  the  positions  of 
actuators  and  sensors  so  that  index  (12)  is  maximized. 

In  general,  the  same  results  can  be  obtained  from  Lagrange’s  dynamic 
equations. 

A  simple  condition  has  been  studied.  Through  analysis,  we  have  derived 
a  set  of  simple  and  practical  formulae  (  7  )  and  (  8  )  that  give  the 
very  positions  at  which  the  effective  damping  force  is  always  dissipating. 
By  applying  formulae  (  7  )  and  (  8  )  to  a  given  missile,  it  is  shown  that 
the  present  position  of  the  elevator  is  suboptimal  and  is  0.06m  away  from 
the  optimal  position. 


NTJMFFICAL  SIMULATION  OP  A  FLEXIBLE  AIRCRAFT 
TAKING- OFF  AND  LANDING  ON  UNEVEN  RUNWAY 


Rui  Yuting 


(Beijing  Institute  of  Aeronautics  and  Astronauti 


PU'T'APY  This  article  does  research  on  the  us-  of  Yen-  equations  to 
set  ur  a  method  of  dynamic  equations  for  flexible  aircraft  sy.  tr  .. 
including  landing  gear  structures.  It  presents  two  types  of  mod -Is 
for  the  calculation  of  runway  surfaces.  After  this,  it  discusses  a 
numerical  simulat^n  method  for  the  attitude,  speed,  acceleration,  and 
dynamic  loading  of  flexible  aircraft  during  the  process  of  taking  off 
and  landing.  It  does  calculations  of  the  effects  on  aircraft  system 
dynamics  of  the  inertial  forces  and  inertial  moments  of  force  of  all 
the  major  components  of  take-off  and  landing  gear.  To  make  a  simple 
example,  this  article  simulates  the  symmetrical  take-off  and  landing 
of  the  KC- 1  35  aircraft. 


I.  Introduction 

When  flexible  aircraft  are  taking  off  and  landing  cn  uneven 
runways,  their  attitude,  speed,  acceleration,  and  dynamic  loading  ar- 
the  main  parameters  for  the  design  of  landing  gear,  calculation  of 
lift,  and  aircraft  capabilities.  Normally,  these  parameters  are 
obtained  through  on-site  measurements  and  high  speed  f notogr ny .  Tn 
costs  are  high,  and  the  time  periods  involved  are  long.  Because  of 
this,  the  use  of  computers  to  do  a  completely  realist!  •  nut  •  ri cai 
simulation  in  order  to  obtain  ti.=s  =  parameters  is  a  pro  tie:.,  worn. 
r~  sear  caning.  ;  1  -3; 

Th“  numerical  simulation  of  ui  rural’*  ta/.e-off.:  an;  landings 
primarily  meets  with  the  following  three  types  of  problems :  1  !  : 

mechanical  model  and  dynamic  equations  for  aircraft  tax— off  an; 
landing  systems;  •'  2  !  a  calculation  model  for  the  surface  of  run*  ays ; 

' 7 '  forms  *0  express  the  effects  cf  t.ae  ground  surface  cn  t ;  r<-  forces 
ani  aerodynamic  f  r  a.  . 
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•  .n  , 


Concerning.  tn®  first  prooxd,  alter  tne  aircraft  lower?  its 
r.g  gear,  the  elastic  fuselage ,  wings,  and  landing  gear  struct . 
a  system  w  i  t  i.  multiple  degrees  of  freedom  and  multiple  otveC*. : 
calculate  the  effects  on  aircraft  system  dynamics  of  the 
ini  forces  and  inertial  moments  of  force  of  all  the  n.  gcr 
r.-  nts  of  landing  gears,  then,  in  that  case,  vie',  we  us-  class, 
nie-.-'l  analysis  methods  t  o  ennur.  e-rate  the  dynamic  equation?  of 
m  ,  we  will  meet  wit:  a  difficult  problem  expressing  complicate 
1  ••?  funkier,-'  wit:,  forms  involving  multiple  numbers  and  scum 
itiss.  At  tr.c  present  time,  in  references  (.1-3/  acing  rcoc.r. 
e  i  y  n  am  i  c  e  1 1  e  c  t of  symmetrical  a.rcratt  t  ari  e— o  i  1  s  an  j  mtiiii 


Method  for  Petting  Up  Dynamic  Equations  in  Aircraft  Systems 


I 

Concerning  systems  compos^ q  of  point  masses  and  holonomic  or 

'  -t 

!  first  order  non-holonomi c  multiple  rigid  bodies,  T.R.  Kane  maintains 

« 

that  it  is  useless  to  solve  dynamic  functions.  They  do  not  have 
f  multiples,  and,  there  is  no  need  to  consider  the  Kane  formula  for 

r-'"Otive  restraints  f4-5l 

F/  +  F*=  0  i*l,  2,  •••,  l—m 

F.  is  the  generalized  motive  force  for  the  system,  /•'•.is  the 
generalized  inertial  force  for  the  system.  /  is  the  generalized 
:•  cord  i  rate  number,  m  is  the  non-holonomic  restraint.  The  method  of 
i  solving  for  F,  and  p»  under  different  circumstances  can  be  seen  in 

references  r 


Pyetem  and  Elastic  Main  Body  1  .  Reference  Body 


J 


bars-  a  1  r--raf*,  whi  1-  taking  off,  landing,  or  taxiing  on  the 
.a;,  jan  \  ~  .-one  i  i-r-d  abstractly  as  main  bodies  and  are  elastic 
:  ,  t  a  1 4  1  :  1  -  body  sys*  »•;:..=• .  Elastic  fuselage  and  wings  are  the 

j  r  1  r.  o  1  t  a  I  r  '  i  1  ■  s  .  be-  landing  gear  components,  ailerons,  and  sc-  on, 
ar-  rigid  t  ody  aux i 11 i ary  components.  Through  generalized  hinges, 

* : .  -  y  ar*-  ■''orneet-d  to  the  main  body.  One  can  see  the  main  body  as 
t .  ;  r.v  tr.-  combination  of  an  infinite  number  of  point  masses.  The  Kan 
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M* 


* 


a 


•ffla 


1 


equation  car.  be  applied  to  this  type  of  multibody  system. 

As  is  shown  in  Fig.  1,  if  we  select  the  Tisserund  reference 
system  i’  to  be  the  floating  coordinate  system  for  the  main  bodies,  ej,  e',. 
are  three  basic  vectors  for  it.  The  position  vector  for  a 
point  p  on  the  main  body  in  the  inertial  coordinate  system 
is  *(*!»  *1*  *j) 

R  =  ',+  r,+  *  (2) 

In  this  equation,  r,  is  the  position  vector  for  the  instantaneous 
center  cf  mass  of  the  main  body.  is  the  position  vector  for  a 

point  corresponding  to  a  point  i  on  the  reference  body,  p  is  the 
elastic  displacement  of  point  P  . 

The  velocity  and  acceleration  for  point  p  are 


R  —  r,  +iX  ('.  +  *)  +  * 

•  *• 

S-h+ixf+ix(»xr)+2oxr+: 


“  is  the  angular  velocity  of  the  reference  body.  r  —  r.+  r,,  t 

•• 

.  :  *  rej  resent  respectively  the  first  ana  second  order  derivatives 
r  ‘v  r-lat  i Ve  coordinate  system  *'  versus  time.  If  we  ns.-u.~f  t hat 


*  h-  r-lat  ive  coordinate  system  e'  versu 
g -  n ■  riu.Zr]  coordinate  of  trie  system  is 

■  v- loc i t ies  set  uj  for  the  system  define 


<7i  *  9*.  "  *  9». 


,  ana 


f.-  2  r«, 

1*1 

l  ■  « 

i-  2 


w: 


r 


modulus 


matrix . 


e±  em^r.  t  c 


(■„  ■„  —  *i)T—  A(q„  Qt*  ••••«)+■«» 

In  this  equation,  i,  is  the  elastic  modulus  for  the  nr.  a i  r  tody  v  iv- 
at  that  in  a  convenient  way,  the  subscripts  •  f'rorr.  7  to  ",  are 
i  r.  such  a  way  that  !?•  corresponds  to  the  f:r^~  cri- r  el  at  1 1  v 
modulus'.  A  i s  an  /  *  /  order  non-singular  matrix.  The  element, 
this  and  the  »  order  line  matrix  A,  are  functions  of  tne  general:, 
coordinate  versus  time.  From  equations  and  viy,  we  :  the 
partial  velocities  for  point  P 


i  = 1 ,  2,  3 
1=4,  5  ,  6 
i  =  7  ,  8  ,  " 

»-«  +  !,  «  +  2 , 


From  the  equation  C^?)“  -  J"K  v.. -5  it  is  possible  to  derive 

the  contribution  made  by  the  elastic  main  body  to  the  generalized 

inertial  force  as  beins 


-  S  )if  +  /T 


i  -  1 ,  2  ,  3 

i  •  4  ,  5  ,  6 

i  “  7  ,  8  ,  •••» 

i  “  *  +  1 »  •••»  / 


In  this  equation,  m  if  the  mass  of  the  main  body,  f  is  t.r 
tensor  of  the  main  body.  ®h  is  th-  g-r.-ral  i  i  mas.  of  : 
The  non-linear  quantity  f*  includes  quantities  s  :■ :  r.d  :  re 
higher  for  the  elastic  displacement  «  and  tr.e  angular  v< 
of  the  main  body.  When  we  are  talking  abou  +  ai rcraf t  in  " 
taking  off  and  landing,  it  is  possible  to  know  the  sub  qua; 

and  *  .  If  we  ignore  the  higher  order  quantities,  then, 
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V  *1  u  r 


x  > 


i  » 


From  equation  )v  v"  “r  '  "  ’ 

possible  to  derive  the  internal  force  contribution  a 
body  to  the  generalized  motive  force,  which  is 


t  h  t 


j  C  lli  'h. 


( F ,) 


j  -mfijq, 

l  0 


i  *  7  ,  8  ,  n 

remainder 


(i 


In  this  equation,  is  thf  *  th  order  generalized  frequency  for  the 
elastic  main  bod^ . 

In  this  system,  the  contribution  of  all  rigid  body  parts  to  the 
generalized  inertial  force  can  be  expressed  as 


i 

(F*).--  ^  *./“/-*.  i  “  1  ,  2,  •••,  | 

i  -  l 


(1C 


In  this  equation,  _  and  *i  are  g-neralized  coordinates  as  a 

*11 

function  of  the  determined  velocities  and  time. 

Or.  the  basis  of  the  basic  formula  (F,)  —  v.rF  i  —  1 ,  2,  — ,  /, 

we  solve  for  all  external  forces  acting  on  the  syc Lem  and  for  the 
contributions  of  all  active  restraining  forces  in  the  system  on  the 
generalized  motive  force.  We  take  that  and  equations  (8),  (9).  and 
(10)  together  and  substi+u+e  them  into  Kane  equations.  We  obtain  the 
dynamics  equation  for  the  elastic  body,  multi-body  system  which  is  the 
m  a ; n  bod y 


,  u,y-H  '  1 

In  this  equat  '  M  |  x  1  ,y  .m-trical  square  matrix.  H  ...  I 

order  line  matrix.  Equation  (6)  and  equation  (11)  form  •  2I  * 

order  non-linear  differential  equation  system. 


ITT.  A  Calculation  Equation  for  Runw  y  Garf aces 

Case  1  i  i  n  the  cdSc  o  1  straight  line  iuiCe— oils  of  £ii  rcr;if  L  f 
landings  and  6rouni3  surface  tauc-off  runs,  it  is  possible  to  reco^ni 
that  the  tires  on  the  three  landing  gears  of  the  aircraft  will  make 
three  mutually  parallel  movements  in  vertical  planes.  Therefore,  wh 
one  is  simulating  the  straight,  line  take-off  and  landing  of  eircra'’ 
it  would  be  a  good  idea  t;  measure,  on  the  surface  of  th-  runway  3 

grid  point  locations  of  heig;4  *•<  1  “  *  •  2,  3»  /  “  1 .  2,  •••,  m. 

If  v;e  look  st  figure  2(a)  ,  2a  in  i F.g  ip  nasrOv  ontigl 

to  the  distance  between  the  two  primary  lauding  gear  on  the  left  and 

right  of  the  aircraft.  The  lengtr.  of  &  should  be  fixed  with  a  view 
to  tr.e  wave  lengtn  of  tne  ups  and  downs  v.  1  tnc  runway.  In  reference 
[11,  the  selection  is  .b  *6icm.  ■  We  set  up  three  one  dimensional 
list  functions  */»  *1/1  xi>  ;*l,  2,  •  ,  m. 

We  use  one  dimensional,  three  poins  equiuisiant  .goiation  values 

to  solve  for  the  h.igr's  for  points  at  die4  —  a  j  x  >,  that  is, 

h,(  3 r  )  t  =  1  ,  2,  3. 
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-P  * 

'•  » 


deflection  forces  whos-  moments  have  effects  on  air  or 


i  n  o  r  d  r  r  x  c 


determine  the  height  of  any  given  location  h\*t  >  )*  ,  ...  >• .  ■ 

the  h  e  i  gh  t  of  a  1 5  x  w  xric  point  i  ■  1  i  2 ,  •••»  5i  i  “  1  ♦  2, 

As  is  shown  in  Fig.  21b),  we  set  up  a 

two-dimensional  equidistant  list  function  and  use  a  two-dimension., 
equidistant  parabola-parabola  interpolation  to  calcuxatt  fc(*.  > ). 


IV.  Ground  Surface  Forces  and  Aerodynamic  Forces 

Forces  having  effects  on  the  tires  can  be  analyzed  into 
supporting  forces  along  the  vertical  direction  Fr,  ,  friction  forces  f , 
and  lateral  fore  es  Fu  ,  which  act  on  tires  when  aircraft  are 
sliding  through  turns.  Ft  is  determined  through  static  pressure 
characteristic  curves  for  the  tires.  Fh  is  determined  from  the 
integral  motion  equations  for  the  rotation  of  the  wheels  of  the  air¬ 
craft  or  from  storage  friction  forces.  Below,  we  discuss  the  solution 
for  the  lateral  forces  Fs 

As  is  shown  in  Fig.  3,  if  we  assume  that  the  direction  of  the 
aircraft's  forward  landing  gear  wheels  cannot  be  controlled,  the 
lateral  force  which  affects  the  main  wheel  on  the  left  and  the  lateral 
force  which  affects  the  main  wheel  on  the  right  are  equivalent.  The 
instantaneous  rate  of  curvature  radius  for  the  track  of  the  movement 
at  a  point  \C'  between  the  left  main  wheel  and  the  right  main  wheel  is 

.  o(*V+*V) 

P  Vw  -V.  .1 


Fig.  3  Force  Analysis  for  Aircraft  Slide  in  Turns  1)  Fron 


3)  Right 
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In  this  equation  Vw  (left)  an’  V .  (right)  represent  respectively 

the  movement  velocities  of  the  left  and  righ+  main  wheels.  Speaking 

in  general  +  e  v  us ,  the  distance  from  the  «.•=  :.-r  of  mass  .C  to  the 

point  C'  must  be  much  smaller  tvnn  P-  .  For  example,  in  the 

BOZI-w  707  aircraft  d-  1.32m,  in  the  POEING  777-94  5  733-94  aircrai  •  ,  d  =1.27 

Therefore,  we  can  know  that  the  r  .  ,.e  of  curvature 

radius  for  the  center  of  mass  ,q  is  equal  to  the  rate  of  curvature 

radius  for  point  C' >•  In  a  natural  coordinate  system,  the  normal 

dynamics  equation  for  the  aircraft  is 

m  ~c-  -  2  F,+  Q-Fh  sin  a 


In  this  equation, < w* •  is  the  mass  of  the  aircraft.  Vc  is  the  velocity 
of  the  center  of  mass  of  the  aircraft-  Q  is  the  horizontal 
component  of  the  main  aerodynamic  force  vector  having  an  effect  on  the 
aircraft.  From  this  equation,  we  can  solve  for 

-Q+F*  sin.)  (1 


as  an  approximation  calculation.  In  the  interval  of  integration  ,Cb, 
we  select  Vclt^'cln,  and  substitute  into  equation 

(13)  to  get  Fs\txF ^,ie 

Aircraft,  during  the  processes  of  take-off  and  landing,  exhibit 
very  large  changes  in  velocity  and  are  affected  by  ground  surface 
forces  which  change  from  time  to  time.  Because  of  this,  aerodynamic 
lift  forces  should  be  calculated  on  the  basis  of  abnormal  aerodynamic 
forces  when  aircraft  are  making  any  movements.  However,  the  problems 
involved  in  carrying  out  calculations  for  non-steady  state  aerodynamic 
forces  when  aircraft  are  making  arbitrary  movements  are  extremely 
complicated.  This  article,  on  the  basis  of  references  (3,7), 
introduces  calculations  of  aerodynamic  forces  from  steady  state  flow 
and  laminar  strip  theory. 


V.  Simulation  Methods 

Let  us  consider  trie  aircraft  system  composed  of  the  elastic 
fuselage,  wings,  forward  landing  gear  and  two  main  landing  gear. 

4  Assume  that  the  wheelbarrow  type  landing  gear  is  composed  of  the  main 

*  column  ,  the  plunger  Bt,  ,  and  the  wheelbarrow  or  trolley  frame  B 

(three  rigid  bodies).  See  Fig.  4(a)*  Assume  that  the  rocker  arm  typ 
landing  gear  with  vibration  reduction  devices  on  the  outside  of  the 
main  column  is  composed  of  the  main  column  •  ,  the  rocker  arm  Bt, 

the  exterior  vibration  reduction  device  tube  Bs  and  the  plunger  Bv 
See  Fig.  4(b).  Let  us  further  assume  that  the  rocker  arm  type 
landing  gear  with  vibration  reduction  devices  inside  the  main  column 
is  composed  of  the  main  column  B i  ,  the  rocker  arm  B.  and  the 
plunger  b5  •  See  Fig.  4(c).  Aircraft  on  which  the  forward  landing 
gear  are  rocker  arm  type  and  the  main  landing  gear  are  wheelbarrow  or 
trolley  type  or  other  types  of  aircraft  are  composed  of  different 
combinations  of  the  three  basic  types  of  landing  gear  described  above 


Fig.  4  Models  of  Landing  Gear 

The  inertial  coordinate  system  <_  (e„  it,  «,)  ps  fiXei  on 
the  surface  of  the  ground.  Its  three  basic  axs-s  run  respective.lv 
elor.g  the  horizontal  direction  of  the  runway,  the  horizontal 
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direction  perpendicular  to  the  runway,  and  the  vertical  direction.  We 
selec^  the  Tisserand  axis  system  «5*  **)  to  be  the 

reference  axis  system  for  measuring  elastic  deformations  in  aircraft. 
Its  three  basic  axes  correspond  respectively  to  the  two  perpendicular 
horizontal  axes  and  the  vertical  axis  of  the  aircraft  reference, 
body.  Let  us  assume  that  x,  2  '  s  the  origin  peir.;  of 

coordinate  system  .  In  coordinate  *  ,  \|i,  6,  'P  ar-.  the 

tnree  Klaylov  a es .  07*  0i*  •■■»  0«  are  the  Generalized 
coordinates  f  "  ‘ir  symmetrical  and  asymmetrical  modalities  of 
aircraft.  si  is  the  travel  path  of  the  vibration  reduction  device  on 
the  forward  landing  gear.  If  one  is  dealing  witn  a  wheexte.r  rev.  cr 
trolley  type  landing  gear,  then,  the  angle  of  rotation  of  iv-_- 

tc  B i  is  p1#  .  Corresponding  to  what  was  said  above,  the  ieft  icon 
gear  quantities  *i  and  0».  .  The  right  main  gear  quant i  '•Ko  are  *s 

end  p,.  .  With  selected  velocities, 

»,  1 

u,  1 

1 

:  -5. 

;  =  ca 

c,c. 


In  this  equation,  5,= sin  o  ,  (?,«cos&,  5,-sin«r,  C.-cost, 

Obviously,  equations  (5)  and  (6)  are  satisfied.  W  —  n  ®  “V— 
in  the  equation  above,  the  square  matrix  /  on  ri0:.‘. 

side  can  be  reversed. 
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On  the  basis  of  the  previous  discussion,  we  solve  for  tut 
contributions?  of  •-•11  the  rigid  body  components  of  the  front  1-  tiding 
gear  end  the  two  mein  lending  gear  to  the  generalized  inert! -si  force . 
We  solve  for  the  contributions  to  the  generalised  motive  force  by 
aerodynamic  forces,  the  ground  surface  forces  affecting,  the  tires,  tic- 
gas  pressure  forces  inside  the  vibration  reduction  devices,  hydraul  ie 
drag  forces  and  storage  friction  forces  as  well  as  trolley  gear  araL 
corresponding  to  the  drag  forces  of  train  column  rotation.  If  we 
assume  that  the  elastic  deformation  of  the  aircraft  ana  the  angu.ar 
velocity  of  the  aircraft  are  small  quantities,  we  ignore  their  small 
second  order  components.  We  consider  the  structural  drag  of  the 
fuselage  and  the  wings,  and,  in  equation  (11),  the  right  side  drag 
quantities  are  ( 0 •  0 »  0,  0,  0,  0,  A,«„  A,«„  •••,  A.«„  o,  •••,  o  )T, 
finally ,  we  get  the  aircraft  system  dynamics  epohicr. 

Af(u,,  U„  . , 


from  equation  v  1  4 ;  ,  we  g- 


( y*  . *  3j)  *  A  (»,.  . •/*) 

Toe  equation  set  of  (15)  and  (16(  is  simultaneously  solved  by 
integration  using  a  computer. 

Inside  vibration  reduction  devices,  there  are  storage  friction 
forces  F,,  <  —  i ,  2,  3.  .  The  subscripts  1,2,5  respectively 

represent  tne  front  landing  gear,  the  left  main  landing  gear,  and  th- 
right  main  landing  gear.  Within  the  area  of  integration  Cb» 


i  -  1  ,  2  .  3 


thfcn,  the  vibration  reduction  devices  are  locked  in.  At  this  tim- ,  w< 
add  he  restraining  conditions 


9f  "•♦iM  " 


We  eliminate  the  "  +2'  *  th  equation  from  the  equation  set  of  (15; 

and  ( 1 6 ) . 


VI.  Sample  Calculation 

This  article  makes  use  of  ATLSP  (Airplane  Take-off  and  Landing 
Simulation  Program)  simulation  programing  to  make  calculations  of  the 
process  of  symmetrical  landings  by  KC-1 35  aircraft.  Initially, 
aircraft  data  was  taken  from  reference  [2].  Runway  data  was  taken 
from  reference  [1],  All  missing  parameters  such  as  landing  gear 
component  mass  and  moment  of  inertia  are  supplied  as  assumed  values. 
The  increment  of  integration  is  taken  as  A*  *  0.001  s ,  .  The 

reliability  of  the  ATLS?  programing  is  supported  by  large  amounts  of 
empirical  measurements  and  experimentation  carried  out  with 
domestically  produced  aircraft.  We  carried  out  a  comparison  of  these 
with  the  calculation  results  from  the  simulation  in  order  to  verify 
it . 


Fig.  5  Aircraft  Angle  of  Pitch  1.  Stationary  Aircraft  Configuration 
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Fig.  6  Vertical  Forces  Influencing  the  Front  W.neexs  ana  gain  m.rcx 
1 .  Stationary  Aircraft  Load  2.  Front  Wheel  5 .  Fain  Wheels 

The  aircraft  assumes  a  landing  speed  of  71.6  m/s .  It  lands  4 
descent  speed  of  1.8  m/s .  3  seconds  after  main  wheel  touchdown, 

flight  personnel  pull  the  stick  5  degrees-  Fig.  5  is  the  aircraft 
pitch  angle-time  curve.  From  this  Fig  we  can  see  that,  5  seconds 
after  the  aircraft  lands,  it  reaches  the  pitch  angle  for  a  station-, 
aircraft.  After  this,  it  oscillates  in  a  decaying  sine  wave  form. 
Fig.  6  is  the  load-time  curve  for  the  load  affecting  the  front  whet 
and  main  wheels.  From  Fig.  6  we  can  see  that,  after  the  aircraft 
first  touches  the  ground,  it  bounces  and  leaves  tne  grouru.  Then, 
executes  a  second  landing.  The  vertical  forces  affecting  the  main 
wheels  at  the  moment  of  impact  are  certainly  not  maximum.  Fig.  7  i 


Fig.  7  Vertical  Acceleration  at  Poinx  C 
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NUMERICAL  SIMULATION  OF  A  FLEXIBLE  AIRCRAFT 
TAKING-OFF  AND  LANDING  ON  UNEVEN  RUNWAY 
Rut  Yuting 

( Beijing  Institute  of  Aeronautics  and  Astronautics ) 

Abstract 

The  method  of  using  Kane’s  equation  to  develop  the  dynamic  equati 
ons  of  a  flexible  aircraft  system  including  landing  gear  mechanisms  is  in 
vestigated  in  this  paper.  Two  kinds  of  mathematical  models  of  runway 
surface  are  given.  Then,  the  numerical  simulation  method  is  discussed  for 
simulating  the  attitude,  velocities,  accelerations  and  dynamic  loads  of  a 
flexible  aircraft  during  take-off  and  landing,  taking  into  account  of  the 
effect  of  inertial  forces  and  moments  of  large  landing  gear  parts  on  the 
motion  of  the  system.  As  a  simple  example,  the  symmetric  take-off  and 
landing  of  KC-135  aircraft  on  the  runway  No.  12  is  simulated.  I  he 
numerical  results  are  only  for  reference. 


